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Abstract 

I derive sharp semiclassical asymptotics of / |e/,(x, y, 0)pa;(x, y) dx dy where e/,(x, y, r) 
is the Schwartz kernel of the spectral projector and u>{x, y) is singular as x = y. I 
also consider asymptotics of more general expressions. 

Introduction 

In the series of papers [IS, Ivr2, Ivr3, Ivr4] devoted to the Sharp Asymptotics of the Ground 
State Energy of Heavy Atoms and Molecules it was needed to calculate Dirac Correction 
Term^^ which in that approximation was equal to 

(0.1) / = = // \e{x,y,T)\^\x-y\-^dxdy 



where e(x,y,r) is the Schwartz kernel of the spectral projector E[r) of the (magnetic) 
Schrodinger operator 

(0.2) ^ = Pjg''{x)Pk - V) . Pj = hDj - f,Vj 

j.k 

r ^ and h — *■ +0 (while either fi —>■ +oo or remains constant). Actually the corresponding 
part of these papers was originally more complicated but it was reduced to the problem 
above. 



Representing Coulomb interaction of electrons with themselves which should not to be counted in the 
energy calculation and should be subtracted from the Thomas-Fermi expression. 
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Then I x h ^ where d is the dimension (c/ = 3 in the above papers) and it was needed 
to prove that 1=1+ 0(/)~^"^+^) with J defined by the same formula but with e(x, y, r) 
replaced by 

(0.3) e;'(x, y, r) {Inhy' I e''''"<^->''«> di 

Jg{y.i)<v{y)+2T 

and with a small exponent 5 > 0; for Magnetic Schrodinger it was needed to prove as < 
only. Expression (0.3) is a Weyl expression for e(x, y, r) for operator with coefficients 
frozen at point y. 

However I believe that the asymptotics of expression (0.1) or more general one is inter- 
esting by itself and that there are a sharp asymptotics. Still my attempts to derive it were 
not very successful and in [Ivr6] I made some claims which I could not sustain at that time. 
So in this paper I just want to bring some degree of the order to this matter. 

I am going to consider a matrix /)-differential operator A[x, hD) and find asymptotics of 

(0.4) JJ ^(-^' ■^)^(-^' ''^)^2(-><')e(y, X, r)^i(y) c/x c/y 

with a matrix-valued function c<j(x,y) such that 



(0.5) uj{x, y) = Q(x, y; X - y) where function Q is smooth in B{0, 1) x B{0, 1) x e(M^ \ 0) 
and homogeneous of degree —n {0 < k < d) with respect to its third argument^-* 

and with smooth cut-off functions ipi,ip2- 

The main part of asymptotics should have a magnitude of /7~'^~'^ and I would like to get 
a remainder estimate 0(h^~'^~'^). 

One can also consider a more general expression 

(0.6) = J J uj{x^ x"')e(x^x^r)V'2(x')e(x^x^r)■■■e(x'",x^r)^^+l(x°)x 

c/x^ • ■ ■ dx"" 

with x™'*"-'- = x-^, ipm+i == i^i etc and 

(0.7) uj{x^, ... , x"') = Q(x\ ... , x'"; {x^-x-'+^}i<^<„) where function Q is smooth in B{0, l)'"x 
6(M'^\0)'"~-'- and homogeneous of degree — (m — 1)k with respect to {x-' — x''}i<,</(<m. More- 
over, 

\D^D:Q\ < Clz^r'^-l^'l-.-lz'"!-^-'^'"' as 5^|zY = l,$^z'^ = 

where x = (x-^, ... , x"'), z = {z^, ... , z'"), etc. 

^) In other words it is Michlin-Caldcron-Zygmund kernel. 
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However I will leave it for another paper since not of all my arguments I was able to 
implement in this case. 

The main part of asymptotics should have a magnitude of l^-d-(m-i)K ^g^g Theorem 1.6) 
and I would like to get a remainder estimate 0[h^~'^~^'^~^^'^). 

I am also leaving for another paper the similar but much more delicate and difficult 
analysis for a 2-dimensional Magnetic Schrodinger operator(0.2) with the trajectories having 
many loops. 

Remark 0.1. (i) To avoid the necessity to cut-off with respect to hD one needs to assume 
that its symbol satisfies 

(0.8) \a{x,i)\-'<C\i\-'^ as|e|>Co 

as a G U''" (one can weaken this condition but I leave it to the reader); 

(ii) One needs to assume that a is semibounded from below which under (0.8) is equivalent 
to 

(0.9) {a{x,C)v, v) > c^Vl^ as |^| > Co; 

otherwise instead of E{t) one should consider E[ri,T2)'= E{t2) — E[ri); I leave it to the 
reader as well. 

This paper consist of two sections: in section 1 I derive asymptotics with the sharp re- 
mainder estimate but with the implicit Tauberian approximation for e(x, y, 0). In section 2 
is I replace it by expression (0.3) without deteriorating remainder estimate for scalar oper- 
ators under mild non- degeneracy condition (theorem 2.19) and for certain matrix operators 
(theorem 2.20(i)) and with some not sharp remainder estimates for other matrix operators 
(theorem 2.20(ii)). I just mention that for larger k we need less restrictive conditions to 
operator. 

1 Estimates 
1.1 Special case 

Let us assume first that u = 1 but relax conditions to ipi: ■■■ •i'm, assuming only that 
ipi, ■■■ ,V^m £ L°°- This is definitely not the case I am interested in but one needs to make 
few clarifications first. Then 

(1.1) L= TrE(r);^2^(r)^3^r)---E(r)^„+i 
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containing m factors E(t). 

Under condition (0.9) it is known (see f.e. [Ivrl]) that if L°° norms and the diameters 
of supports ip, ipi are bounded, then 

(1.2) |||^/'E(r)V^i|||i < Ch-'' as |r| < c 

where f-f^ and |||.|||j^ denote operator and trace norms respectively. Then since an operator 
norm of E{t) does not exceed 1 I conclude that |/m| < ch~'^ . So 

(1.3) If e and is given by (1.1) then |/^| < Ch-'' . 
Further, let us assume that 

(1.4) a(x, ^) is microhyperbolic on energy level 0. 

Then as (1.4) is fulfilled on supports of ip, ipi it is known (see f.e. [Ivrl]) that 

(1.5) \ltP{E{r)- E{T'))^^i\l<C{\r-r'\ + hT-^)h-'' as |r| < ei, |r'| < c. 

Here and for a while 7~ x 1 but I want to keep a track of it. 

Since this property holds under wider assumptions than microhyperbolicity, I will as- 
sume so far only that (1.5) holds. 

Then 

(1.6) I Tr' ((E(r) - E(rO) V2^(r2)^3f (ts) ■ ■ ■ E(r^)7/wi) I 

also does not exceed the right hand expression of (1.5) as |r| < ei and therefore due to the 
standard Tauberian arguments (second part, see f.e. [Ivrl]) the following inequality holds: 

(1.7) |Tr'(^(E(0)-/j-^y"° E,^,-i,(xr(f)^0)c'r)^2f(r2)^3^(r3)---E(r^)^„+i^| < 

where I use my standard notations x X the future and x(t) = x(t/T) etc (see f.e. 
[Brivr]). Here and below Tr' is the "scalar trace" of the operator, and does not include 
taking matrix trace tr. 

Here and below U{t) = e'^ '^^ is the propagator of A and u{x,y, t) is its Schwartz' 
kernel. 

So with 0( 7""^/?-^"'^) error one could replace one copy of E(0) in 1^ by its standard 
implicit Tauberian approximation 

(1.8) h-' f F,^,-.^{xT{t)U{t))dT 

J — oo 

and in by the virtues of the same arguments I can do it with another copy of E(0). Therefore 
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Proposition 1.1. Under conditions {1-5) with an error 0{T~^h^~'') 1^ is equal to 

(1.9) h-^^Tr' [ F,^,-iJxT{ti)U{t,)^2xAt2)Uit2)^P3 - ■ ■ U{t^)^^^^)) dr 

with t = {ti tm), T = (ri Tm). 

Note that here one can take any T E [Ch^^^ , c] (but then an error depends on 7"). 
Further, note that as dist(supp ^/'j, supp > (cq + e)7 where Cq here and below is the 

upper bound of the propagation speed on energy level and x'""^-^ =^ x-*-, expression (1.9) as 
m = 2 or the similar expression as m > 3 become negligible and I arrive to 

Corollary 1.2. If in frames of proposition 1.1 dist(supp ^/'y, supp ^/'j+i) > [cq + €)T for some 
j = 1, ... , m then does not exceed CT~^h^~'^ . 

1.2 Smooth case 

The next step is to assume that is a smooth function. Without any loss of the generality 
one can assume that uj is also compactly supported (since are). Then from 

(1.10) uj{x' x'") = j uj{y^ y'")5(y'-x^ y'"-x'")c/y = 

one arrives to 



;i.ii) L = j y"')^2(y' y'^)dy'---d/ 



with J2{y^, ... ,y'") defined by u; = 1 and ipj{x) redefined as tljj{x)6[y-' — x) where here and 
below 6{x) = 6[xi) ■ ■ ■ 6{xd). Then I immediately arrive to 

Proposition 1.3. Let uo and ipi, ... ,1^^ be smooth functions and let condition {0.9) be 
fulfilled. Then |/^| < Ch^'' . 

Remark 1.4. As m = 2 and uj,ipi,ip2 £ I/2I < Ch~^ obviously (it follows from the 
estimate |||?/'E?/'|||2 < Ch~'^^^ where |||.|||2 is the Hilbert-Schmidt norm). Can one prove the 
similar result for m > 3? 
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Proposition 1.5. Let uj and ipi, ... ,11)^ he smooth functions and let conditions {0.9) and 
{1-5) be fulfilled. Then 

(i) with an error 0{T^^h^~'^) 1^ is equal to 

(1.12) I^ = h~'^ [ I u{x^ x'^)F,^,-.JxT{h)u{x\x\tr)'4,2{x^)xT{t2)^ 

u{x\ x\ t2)^3(x') ■ ■ ■ ty(t^)^^+i(x"'+^)) dr dx'--- dx"' 

with x'"+i = 

(ii) Further, if d\st{supp ipj , supp tpj+i) > (co + e)7" for some j = 1, ... , m then does not 
exceed CT^'h'^^ where so far 7" x 1. 

1.3 Singular homogeneous case 

Theorem 1.6. Let conditions {0.9) and {0.7) be fulfilled. Then < C/)~^"('"~^)'^. 

Proof. Let us replace Q(x, z) by Q(x, z)/3(z-^/7i) ■ ■ ■ /3(z"'/7m) where 7y > h and /3, /3 are 
functions (on R'') similar to x. X respectively. Then similarly to the analysis of the smooth 
case one can estimate the contribution of such partition element to by 

(1.13) C/j-^(7i • • ■7-)~'(7i + ■ ■ ■ + 7-)'"'^ 

and summation with respect to 7^ > 7 = /) results in the value of this expression as 7j = 7 
and the total estimate becomes what is claimed. 

However one needs to consider the other partition elements when some of (3{z-' /^j) are 
replaced by P{z-' /j). So we get "sandwiches" consisting of the factors 

e(x\ x'+\ r)/3(z'^+V7/c+i) ■ ■ - P{z^ /lj)e{xJ , x^+\ r) 

with J > k and in between them factors (3{z'^/'j). 

Let J be the set of indices appearing in /3{z'^/'j) (for a given type of a "sandwich"). One 
can see easily that the contribution of each "sandwich" to does not exceed 

where r is the number of factors of each type. Then after summation with respect to 7^ > 7 
one gets the same expression with 7^ = 7 i.e. C/)"'^''7''("'"''^+('^~'^)(''~^^ = C/)"'^''7"''('"^^)+'^(''~^) 
which is exactly what we want as 7 x /j. □ 
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It immediately follows from the proof a stronger condition 

Proposition 1.7. Let conditions {0.9) and {0.7) be fulfilled. Then replacing Q(x, z) by 
fi(x, z)/9(z^/7) ■ ■ ■ /5(z'"/7) results in the error not exceeding 

(1.14) Ch-''-^'^-^^^j-\ 
Now let assume instead of condition (1.4) or (1.5) that 

(1.15) a(x, ^) is microhyperbolic on energy level and microhyperbolicity directions are (at 
each point) £^ ■ with = i^{x, ^). 

Proposition 1.8. Let conditions {0.9), {0.7) and {1.15) be fulfilled. Then replacing Q.{x, z) 
by f2(x, z)/3(zV7) ■ ■ • /9(z'"/7) results m the error not exceeding 

(1.16) (^^l-d-{m^2)K^~l~n^ 

This is equivalent to taking T x j in {1.8) and plugging Schwartz kernel of it instead of 
e(x,y, 0) in the definition of 1^. 

Proof. Proof follows from the combined arguments of the proofs of Theorem 1.6 and Propo- 
sition 1.1; in this case one needs to consider only "sandwiches" containing at least one factor 
P{x^'/'-fj) with 7j > 7 which accounts for a factor h/'jj and summation with respect to par- 
tition results in an extra factor h/'-f. □ 

So one needs to study expression (1.12) with some T = T*;l remind that the remainder 
estimate contains factor T*~^. One can decompose XT*{t) into the sum of Xf{t) and XT{t) 
with T running between T and 7"* and also one can take T = Ch. Then expression (1.12) 
becomes the sum of the similar expressions with Xt(0 (with T = T*) replaced by 4'jTj{t) 
where either (pj = x and T < Tj < T* or (pj = x and Tj = T. 

In this expression as (pj = x one can replace j_^{...) dr by (...)|,-=o simultaneously 
replacing h~^XT{t) by it~^XT{t) = T~^(f)r{t) with (j){t) = it~^x{t)'^ so we get a modified 
expression (1.12) with r factors Xf{tj) and tj snapped to for j ^ J, r = and integration 
over R~('"^'') and (m — r) factors (pT{tk), k ^ J; furthermore, factor /?"'" is replaced by 

Proposition 1.9. Let conditions {0.9) and{1.15) be fulfilled and letuj be a smooth function, 

(1.17) a;= 0((|x^-x2| + ... + |x'"-x^|)'^). 
Then = 0{h^-'') as K > 1 and 1,^ = 0(/)^"^| \ogh\) as K = 1. 
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Proof. Proof follows from the combined arguments of the proofs of Theorem 1.6 and Propo- 
sition 1.1 like in Proposition 1.8. Here however the main contribution (as K" > 1) is delivered 
by zone {|x^ - x^l H h |x'" - x 1}. □ 

One can consider certain generalizations but I will do it later. 

2 Calculations 

Now our purpose is to go from implicit Tauberian expression (1.12) to more explicit one. 

2.1 Constant Coefficients Case 

Let us first consider case A{x, ^) = A[^). In this case 

(2.1) e(x, y, r) = {Inh)-' J e'''"'<"-^'«>E(0 
where E(^, r) is the matrix projector corresponding to Then 

(2.2) = (27r/?)-'^'" y j u{x^ x'")E(e\0)---E(e'",0)x 

(2.3) From now and until the end of the paper I am assuming that m = 2. 
Without any loss of the generality one can assume that either u;(x, y) is of the form 

(2.4) uj{x,y) = Q{^{x + y),x-y). 

or it is of the same singular type as before but multiplied by {x^ — Yk)- However in the 
latter case (under microhyperbolicity condition) one can apply a Tauberian approximation 
for e(x,y, r) equal with the remainder estimate 0{h^~^\x — y\~^) (in the same trace class 
as before) which leads to / ~ with the sought remainder estimate 0{h^~'^~'^). 
In the former case (2.4) we get 

(2.5) / = l2= [ J{x) dx, 
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where 



(2.6) 



J{x) =2{2nh)-^' / // Q(x,z)E(e,0)E(r/,0)e' 




dzd^dr] = G{x)h 



with 



(2.7) 



G(x) 




77)E(e,0)E(r/,0) d^cyr/, 



and 



(2.8) 



Q(x, C) = 2(27r)-2rf / Q(x, z)e'<"'^> c/z 




One always can take Q having a compact support with respect to x (since we had originally 
cutoffs iJi{x^), ... , V^m(x'"). 

Remark 2.1. (i) One can easily generalize (2.5)-(2.8) to m > 2. 

(ii) Integral (2.8) converges as |z| < 1 since k < d. On the other hand it defines a 
distribution with respect to ( which is positively homogeneous of degree k — d and also is 
smooth as C 7^ 0; thus Q G L]^^ and (2.7) is well-defined. However generalization to m > 2 
is not that easy. 

2.2 General Microhyperbolic Case 

Note first that due to the microhyperbolicity condition (1.15) one should take 7~ x 7 as 
m = Otherwise as 7" G [Ch^~^ , T*], T* is the small constant, the contribution of 

[T/2, T] U [-7,-7/2] would be negligible. 

To calculate u let us apply the successive approximation method on the time interval 
[— 7, 7] with h^~^ < 7. Then plugging the successive approximation into any copy on 
that interval we arrive to an error in u in the trace norm equal to 0[h-'^[T^ / hf) where 
n is the number of the first dropped term (starting from 0). This leads to the error in / 
0(/7~''~'^( 7^//))"7~'^) as 7 > 7. Since under microhyperbolicity assumption (1.15) we need 
to consider only 7x7, the error is 0{h-''{T^/h)"T-''). However if we just take u = 
then we get an error 0{h^^'^ T~^~^). 

Finding 7 from the equation 



-1 



And Tj X Ix-* — x-'+-'-| in the general case. 
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we get 



(2.9) 



T = /,("+i)/(2''+i) 



(which is greater than ^ with 5 > 0) and this leads to an error 

(^2.10) Q^^l-d-(n+l)(K+l)/(2n+l)j 



Proposition 2.2. Let conditions {0.9), (0.7) and {i-15) be fulfilled. Then 

(i) Using successive approximation as \ t\ < T given by [2.2) and taking u = otherwise we 
get I with an error given by [2.10). 

(ii) In particular this is the sharp remainder estimate 0[h^~^~'^) as 



in particular, as n >2 one can skip all perturbation terms and get the same answer [2.4) — 



On the other hand, if we cannot skip some term, then this is given by the same formulae 
(2.4)-(2.7) as before but with the factor instead of h~'^~'^ and with Q replaced by 

Qs positively homogeneous of degree — k + s (provided these formulae have sense!). Then 
as long as s < K one can see that these terms are less than the remainder estimate and we 
arrive to 

Proposition 2.3. Let conditions [0.9), [0.7) and [1.15) be fulfilled. Then 

(i) As K> 1 formulae [2.4) — [2. 7) provide an answer with the remainder estimate 0[h^^^~'^). 

(ii) As K < 1 formulae [2.4) — [2.7) provide an answer with the remainder estimate 
0[h2^^^'^^~'^~'^~^) with arbitrarily small exponent S > 0. 



2.3.1 Assume first that uj = 1 and 1^1,4^2 are smooth functions. Then one can rewrite 
(1.9) with m = 2 



(2.11) 



K>[n + l)/n; 



[2.7). 



2.3 Scalar Case 



Let us completely analyze the case of scalar operator A. 



(2.12) 
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with T = T* which is the largest value for which remainder estimate 0{T''^h^~^) for the 
standard asymptotics was derived; here T* x 1. 

If we replace some copies of Xri^k) by XTk{tk) with Ch < < T* then one can replace 
also operator /j"^ dn by T'^^...) \r,=o and x by it'^X- 

If we do it with both k = 1,2 then we get a term 0{h~'^) (the better estimate is actually 
possible) and the summation with respect to all partitions with respect to T"i, T2 results in 
0{h~'^\ log/?p) which differs from the proper estimate by | \ogh\^ factor. If we replace some 
copies of Xritk) by Xh{tk) then we do not make a transformation with respect to these 
factors but we gain a factor h due to the size of the support. So after summation with 
respect to partition we arrive to estimate 0{h~'^\ \ogh\^"'') for / where r is the number of 
Xh{tk) factors. 

On the other hand expression (2.12) is equal to 

(2.13) h-^Jr / Ft,^,-i^^^t2^,-iJxT{ti)Mt2)i^i'^2.HU{ti + t2)))dT 

with ijjt = U{t)ipU{-t). 

Applying standard approach we arrive to 

(2.14) X~^x„/)"''+" 

n>0 

where X = X2 is defined by (1.12). 

Let us replace in (2.13) ^/'2,ti by ^2- Plugging ti2 = \t ± z, ti^ = t ±t' we arrive to 

(2.15) /)"^Tr / (f pT{t,T)^ii)2U{t)e~''^'''^ dt)dT 
where pjit, r) = p{t/T,T), r < 

(2.16) p{t,T) = -n-H-' f XT{lt + z)xT{l:t-z)z-'sm{h-'TzT)dz 
is Co°([— 2, 2]) and one can prove easily that 

(2.17) \d^{p{t,T)Tf{t/2))\ < CUl + \r\Th-')-'- Vm, n Vr ^ 0. 

Then due to (2.17) only zone {|r| < h^^^} gives a non-negligible contribution to this 
error and due to the microhyperbolicity condition there | Tr?/;i?/'2L/(t)| < Ch~^{l + \t\h^^)~'^ 
which together with (2.17) implies 
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(2.18) Under microhyperbolicity condition (1.4) expression (2.15) is equal modulo 0{h^-'') 
to the same expression with p replaced by x^{t/2). 

On the other hand, if we replace ip2,ti by ip2,ti — = fiV^2,ti ^^en we can apply the same 
transformation as before just getting rid of one factor and integration with respect to 
Ti, which simply snaps to 0, resulting in expression, similar to (2.15) but with ptp2 replaced 
by 

(2.19) p'it.T.x) = {2n)-'i [ XT{z)xT{t-z)e''-'^^ii,dz 

Jr 

which satisfies inequality similar to (2.17) 

(2.20) |9f"p(t, r)| < C„(l + |r| Th-^)-"" Vm, n Vr ^ 0. 
and therefore 

(2.21) Under microhyperbolicity condition (1.4) this new (2.15)-type expression is 0(/)^~'^). 

So, we are left with expression (2.15) with p(t) = x^(t/2) but due to the standard 
theory we get modulo 0{h^~'^) expression 

(2.22) JriPi^P2E{0) = {27rhy [[ ^Pii;2dxdC 

J J{a{x4)<0} 

So, I is given by (2.22) modulo 0{h^^'^) and therefore 

(2.23) xo = (27r)"^/ // ^iip2dxdi in (2.14). 

J JJ{a{x,^)<0} 

2.3.2 Then in the general smooth case we get 

Proposition 2.4. Let uj and ipi, ... .tp^ be smooth functions and let {0.9) and microhy- 
perbolicity condition {1-4) be fulfilled. Then with an error 0{T~^h^~'^) where T x 1 here 
decomposition {2.14) holds with 

(2.24) xo = (27r)~^ /"/" uj{x, x)iPi{x)iP2{x) dx 

J J{3{x.i)<0} 

Proof. Follows from the standard decomposition (l.lO)-(l.ll). □ 
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2.3.3 Consider now the case of singular homogeneous u. First, let us consider Xy defined 
by (1.12) with u = 1 and ipi,ip2 replaced by ^/'i.-y, ^/'2,7 which are some smooth functions 
scaled at some point z with the scaling parameter 7 G {h^~^ , h^). To have microhyperbolicity 
condition sustain scaling we replace it by (1.15). Then (2.14) implies 

(2.25) ^' ^ >^nmh~''+"l"'~"^'' 



n,m>0 



and obviously 



(2.26) {271)-' [ [ ^i,,(x)^2,,(x) dx c/e ~ V x: 

JJ{a{x.O<0} ^>o 



One can see easily that in (2.25) terms with m = would be the same for operator = 
3o{z, hD) where ao(x, ^) is the principal symbol of A; this z is not necessarily the original 
one, but distance between them should not exceed C7; similarly in (2.26) term with m = 
coincides with the left-hand expression with a(x,^) replaced by a(z,^). 

What is more, under condition (1.15) integration with respect to x is not needed, so all 
these results would hold (without factor in the decomposition and estimates) without 
it; thus one can take z = x (or y, does not matter). 

Thus we arrive to 

Proposition 2.5. Let X' he defined by {1.12) with u = 1 and ipi, 1^2 replaced by ipi^^, ip2,-y 
which are the same smooth functions scaled at some point z with the parameter 7 G 
{h^~^,h^). Let X°' be defined the same way but with U{t) replaced by ty°(t) = e'^ 
where A° = a{z, hD) and later z is set to x. Then T - 1°^ = 0{h^~'^-f'') 

Now we can calculate X in the scalar case: 

Proposition 2.6. In frames of proposition 2.5 as uj satisfies {1.7) and k > X — X° = 
0{h^-'-'^) where X^ is defined for constant- coefficient operator obtained by freezing coeffi- 
cients of A at point x (or y, does not matter). 

Proof. Consider three zones: {|x — y| > 71} with 71 x ^ {7 < |x — y| < 71} with 70 x h^~^ 
and {|x — y| < 7}; then the contribution of the first zone to the reminder for X and X° does 
not exceed C/?^-''7f ^"'^ = 0(/)^"^~^) (while main parts are 0); in virtue of proposition 2.5 
and decomposition of subsection 1.2 the contribution of the second zone to X — X° does not 
exceed 0{h^-''-i~'') = o(/?i-^-^). 

In the third zone one can apply the method of successive approximations resulting in 



m+n+k>l 
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However since the final answer does not depend on 7 only terms with 2m = n are posed to 
survive just resulting in (x+ o[l)^h~'^~^^~^. □ 

Summarizing results of section 1, proposition 2.6 and formulae (2.5)-(2.8) we arrive to 

Theorem 2.7. Let A be a scalar operator satisfying conditions {0.9) and {1.15). Then 

(2.27) / = j J{x)Mx)Mx) dx + 0{h^-'~^), 
where 

(2.28) J{x) =2(27r/))"2d ffj ^^^^ q^q^^^ ^^^^^^ Q)g//,-i(z,5-„) ^^^^^^ = 



J J J {a{x,i)<Q, a(x,r))<0} 



with 



(2,29) G(x)^//E(x,COMx,?-,)E(x,,,0).e., 



'{a(x,^)<0, a(x,»?)<0} 

and Q zs defined by {2.8). 

Remark 2.8. (i) Alternatively one can prove this theorem using oscillatory integral repre- 
sentation of u{x, y , t) as \ t\ < T = e. 

(ii) Alternatively one can replace one or both copies of x in E{x, ., .) or in a{x, 0) by y. 

Definition 2.9. We will refer to formulae (2.27)-(2.29),(2.8) as to standard Weyl expression 
even in the matrix case. However in this case the third parts of (2. 27), (2. 28) should be 
skipped. 

2.4 Schrodinger operator 

Now my goal is to weaken and eventually to get rid off microhyperbolicity condition for 
scalar operators. I start from the Schrodinger operator. 

For a Schrodinger operator condition of microhyperbolicity (1.15) means that 

(2.30) V > eo. 
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If this condition is violated let us introduce scaling functions p{x), 7(x) in the usual way 
7 = e| \/| and p = 7"^/^. 

Then, the contribution of 6(x,7(x))^ to the remainder does not exceed 

(2.31) C(/)/P7)^"'^"''7~'' ^ C/)^"^~>'^~^"''7^~^ 
with p = p[x) and 7 = 7(x) and then the contribution of zone 

(2.32) {(x,y) : |x-y| <e7(x)} 

(where automatically 7(x) x 7(y)) to the remainder does not exceed 

(2.33) Ch^-''-'' J p'^-^+'^7-^ dx 

and with p = 7"^/^ here it becomes 

(2.34) C/?^^'^-" J 7(^-3+K)/2 

obviously, it is 0{h^~^~'^) provided either d + k, >3 or 

(2.35) |\/| + IVUI > eo. 
and d + K > 1 (which is surely the case). 

Remark 2.10. (i) Note that (2.35) is microhyperbolicity condition (1.4). 

(ii) Actually one should take p7 > Ch and thus to add Ch^^^ and Ch^^^ to p,7 respectively 
(but it does not affect our conclusion due to the standard fact that if p7 x /? then h^ff x 1 
and condition (2.35) is not needed. 

Consider now the complement of zone (2.32). Let us redefine there 7(x) as 7(x,y) = 
||x — y| and in this zone condition (2.35) is not needed as one can see easily after rescaling 
6(x,7(x,y)) to 6(0, 1) due to proposition 1.5. 

Therefore as 7 > 7(x) the contribution of 6(x,7)^ \ {zone (2.32)} to the remainder 
does not exceed the same expression (2.31) with p = 'j'^^'^. Then the contribution of the 
complement of zone (2.32) to the remainder does not exceed 

(2.36) C/)^-^-" [[ |x - y |(^-i+-)/2-i-rf dx dy. 

J J {|x-y|>e max(7(x),7(y))} 

One can see easily that expression (2.36) is 0{h^~'^~'^) as c/ + k > 3 (so this case is already 
covered) . 
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Further, expression (2.36) does not exceed expression (2.34) with 7 = 7(x) and expres- 
sion 

(2.37) Ch'-'-'' J{\\og^{x)\ + l)dx 

as d + K < 3 and d + k = 3 respectively and both these expressions are 0{h^-''-'^) under 
condition (2.35). 

Again we get 0( provided either d + k > 3 or condition (2.35) is fulfilled. So, 
we arrive to 

Proposition 2.11. Consider Schrddinger operator. Let either d + n > 3 or condition 
{2.35) be fulfilled. Then the standard Weyl asymptotics {2.27} - {2.29), {2.8) holds with the 
remainder estimate 0{h^~^~'^). 

This completely covers the case d > 3. Furthermore, after proposition 2.11 is proven, 
we can introduce scaling functions 7 = p = e(| \/| + iVl/p)"^/^ + Ch^^^ and then applying 
the same arguments we arrive to 

Proposition 2.12. Consider Schrddinger operator. Let either d + k > 2 or condition 

(2.38) IVI + |V\/| + |VV| > e 

be fulfilled. Then the standard Weyl asymptotics {2.27) — {2.29), {2.8) holds with the re- 
mainder estimate 0{h^~'^~'^). 

This completely covers the case d = 2. As c/ = 1 we get the required remainder estimate 
under condition (2.38). 

Now, combining this with the arguments of the proof of Theorem 4.4.9 of [Ivrl] we get^-' 

Proposition 2.13. Consider Schrddinger operator with d = 1, n > 0. Then the standard 
Weyl asymptotics {2.27) — {2.29), {2.8) holds with the remainder estimate 0{h^~'^~'^). 

Remark 2.14. Actually all above results hold as /t = as well with the singular exception 
of d = 1 when the remainder estimate 0(1) is recovered under condition 

(2.39) El^x\/I>e; 

without it the remainder estimate is 0{h'~^) with arbitrarily small 6 > 0. 
I am leaving easy details to the reader; see also the proof of Theorem 2.19. 
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2.5 Scalar Case. II 

2.5.1 Let us consider general scalar operators. 

Remark 2.15. (i) Actually instead of condition (0.9) one can make a cut-off with respect to 
replacing functions ipjix) by pseudo-differential operators tpj{x, hD) with smooth compactly 
supported symbols; 

(ii) Alternatively we can replace £^(0) by E[t,t') = E(t) — E[t') with conditions satisfied 
for a — T and a — r' instead of a. 

(iii) Alternatively we can replace £^(0) by 

(2.40) E'(r) = / E{0,T')^{T')dT' 

Jr 

with smooth function (p s.t. J^ip^r') dr' = 1. 

In all these cases obvious modifications of the final formulae are needed. 

Now we can introduce scaling functions 

(2.41) 7(x, = e(| V^a|2 + |a|) + Ch'/\ p{x, = ^'/'{x, 

and repeat arguments of the previous subsection; then expression (2.33) will be replaced 
by Ch^-''-^M with 

(2.42) M = J dxd^ X J (iV^al^ + |a|)'""'^/' dxd^ 

(in zone {p7 > Ch}). Therefore we arrive to the remainder estimate 0{h^~'^~'^) provided 
M = 0(1) as now integral in M is taken over 6(0, 1). 

This is definitely the CclSG clS Kj > 3. Assume now that microhyperbolicity condition 
(1.4) is fulfilled. Then M = 0(1) as k > 1; otherwise this condition becomes 

(2.43) J iV^al""^ d/i < oo as < k < 1, j | log | V^a] | c/// < oo 

with E = {a(x, ^) = 0} and dp = dxd^ : da measure on Z. 

Thus we arrive to the following generalization of proposition 2.11: 

Proposition 2.16. Let A be a scalar operator satisfying condition {0.9). Assume that the 
uniform version of condition^^ 

(2.44) r a = V^a = =^ rank Hess^^ a > r 



I.e. \a\ + |V^a| < e implies that Hess^j a has r eigenvalues which absolute values are greater than e. 
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is fulfilled. Then 

(i) As r+hi > 3 the standard Weyl asymptotics {2.27) — [2.29), [2.8) holds with the remainder 
estimate 0{h^~''^''); 

(ii) Under condition {1.4) r + n > 1 the standard Weyl asymptotics {2.27) — {2.29), {2.8) 
holds with the remainder estimate 0{h^~'^~^). 

Proof. In contrast to standard asymptotics we need to consider not points (x,^) but pairs 
(x, ^',y,ri) and the pure standard arguments work in zones 

(2.45) {(x, i; y, r,):\x-y\< e7(x, 0, 1^ - < ep{x, 0} 

where also 7(y,?7) x l{>(,C) ^"^^ p{y'V) ^ P{^'0- Analysis in the complimentary zone I 
postpone until the proof of theorem 2.19 where it will be done in more general settings. □ 

Now introducing scaling functions 

(2.46) 7(x, = e(l V^.^ap + |a|)^/' + p(x, = l{x, 

and repeating the same arguments we arrive to the following generalization of proposition 
2.12: 

Proposition 2.17. Let A be a scalar operator satisfying condition {0.9). Assume that the 
uniform version of condition (2.44)^ is fulfilled. Then as r + k > 2 the standard Weyl 
asymptotics {2.27) — {2.29), {2.8) holds with the remainder estimate 0(/)"^~'^~'^). 

Again, combining this with the arguments of Theorem 4.4.9 of [Ivrl] we arrive to the 
following generalization of proposition 2.12 

Proposition 2.18. Let A be a scalar operator satisfying conditions {0.9) and {2.44)i (^n-d let 
K > 0. Then the standard Weyl asymptotics {2.27) — {2.29), {2.8) holds with the remainder 
estimate 0{h^~''~'^). 

2.5.2 Now we can prove our main result for scalar operators: 
Theorem 2.19. Consider scalar operator. Let conditions {0.9) and 

(2.47) „ Yl I^M ^ ^0 

0<k<n 

with some n be fulfilled. Let u satisfy {1.7) and n> 0. Then the standard Weyl asymptotics 
{2.27) - {2.29), {2.8) holds with the remainder estimate 0{h^~''~'^). 
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Proof Part L In this part of the proof we consider at each step only zone (2.45) where 7 will 
be defined in different ways later. Treatment of the complementary zone will be described 
in Part II. 

So, we proved the statement of the theorem under condition (2.44)i which is equivalent 
to (2.47)2. 

Let us apply induction with respect to n. Assume that under condition (2.47)n required 
estimate is proven. 

In the general case (without condition (2.47)n) we can introduce scaling functions in the 
manner similar to (2.41): 

(2.48) „ i{x,0 = e{ J2 |V5^a|'^/("-''+^))^''+''/'' + CM"+^)/("+2), p{x, = 

0<k<n 

with N = (n + 1)!. 

Therefore under assumption of induction we get again remainder estimate Ch^~^^'^M 
with M given by (2.42) where this time the right-hand expression becomes 

(2.49) „ M = j ^{''-"-2)/("+i) dxdi; 
under condition (1.4) this expression becomes 

(2.50) „ Mx /" 7(-"^)/(''+^) X f (V |Vja|i/(''-''+^))"~'c//i 

■^^ l<k<n 

which is 0(1) under assumption |V^^"^a| > eo (as lower order derivatives with respect to ^ 
are close to 0). This is exactly condition (2.47)n+i. 

So, now we have a proper estimate under condition [2.47)n+i instead of {2.4'T)n but now 
we also need condition {i-4). 

Without condition (1.4) we would need something different; f.e. ignoring integration 
with respect to x one should assume that rank(V"''"^a) + k > n + 2 where the rank of multi- 
linear symmetric m-form G is d — dim Ker G; Ker G = {x : G{x, x^, ... , x"^) = Vx^, ... , x'"}. 
This is rather unusable. 

Instead I want to weaken condition (1.4), replacing it by 

(2.51) „+i,„ I'^'Ml > 

2<j<n+l, l:m+j:{n+l)<l 

for some m > which is not necessarily an integer. Obviously in our assumptions (1.4) 
coincides with (2.51)n+i,i. 
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Let us run a kind of nested induction. So, let us assume that under conditions (2.47) n+i 
and (2.51)n+i,m remainder estimate 0{h^^^^^) is proven. 
Now we can go to something similar (2.46): 

(2.52)„,„ 7(x,0 = e( E |Vjv:ar^-)'/'' + 7, ^ = cM''+^)/('"+''+2)_ 

k,l:k:n+l:m<l 

n + 1 



p(x,0 = 7^'"™+'n^.O. % 



(m + 1)(a7 + 1) - (m + l)/c - (n + !)/■ 
Then we recover remainder estimate Ch^^'^^^M with M defined by (2.42) which is now 

(2.53) „„ M X y dxd^ X y p-{"+i)/{-+i)+(-i) dxdC 
Under condition (2.47)n+i we can assume without any loss of the generality that 

(2.54) a{x,0= E bjix.aCi"-'-'. bo = 1, b, = 0; 

0<j<n+l 

we can always reach it by change of coordinates and multiplication of A by an appropriate 
positive pseudo-differential factor. Then 

(2.55) px l^il+p, p = 7(x,0('"+^)/("+^), 7= |Vj,V>/(''+^)'+7. 

J,k,l:{k+J):n+{l:m)<l 

Then 

(2.56) „„ M X y dxde X J dxd^' 

(with an extra logarithmic factor as the power is 0). Then M = 0(1) as 

(2.57) (m+ l)K/(n+ 1) > 1. 

Moreover, M = 0(1) provided there exists {J, k, I) with \S/^,S/'^bj\ > eo and either k > 1, 
{k + j - 1) : n + I : m < 1, Sk+j-ij < 1 or I > 1, {k + j) : n + {I - 1) : m < 1, Sk+jj-i < 1. 
Therefore one can derive easily 

(2.58) If remainder estimate 0{h^~'^^'^) holds under condition (2.51)n+i,m' ioT every m' < m, 
then it also holds under condition (2.51)n+i,m- 



2 CALCULATIONS 



21 



On the other hand there exists a discrete set {m^}u^i^2,... with mi < m2 < ... such that 
if condition (2.51)n+i,m is fulfilled for m = m,, then it is fulfilled for all m G {m^, m^^i) as 
well. 

This justifies induction with respect to m running this set and therefore remainder 
estimate 0{h^~'^~'^) holds under condition (2.51)n+i,m no matter how large m is. However 
if m is large enough, condition (2.57) is fulfilled and we do not need condition (2.51)n+i,m 
anymore. 

This concludes induction with respect to n. □ 

Proof Part IL However in contrast to standard asymptotics we need to consider not points 
(x, ^) but pairs (x, C,', y, rj) and the pure standard arguments work in zone (2.45). 

It follows from the standard theory that if Qx and Qy have symbols supported in 
e(Px, 7x)- and e{py, 7y)- vicinities of (x, ^) and (y, rj) respectively then 

(2.59) WQ.EQyWl < Ch-'{p,^,Y/\py^yY/^ 

and moreover, if either |x — y| > eo7x or |^ — "/yl > eoPx then 

(2.60) WQ.EQyWl < Ch^~'{p,^;)"^-\py^yY'\ 

Surely the same will be true with (x, ^) and (y, 77) permuted. 

Then contribution of such pair to the error estimate does not exceed 

(2.61) Ch'-'{p.^.Y/^-\py^yY''\x - yr 
if k - y| > eo7x- 

Otherwise contribution of the pair ipxQx and Qy to the error estimate does not exceed 
Ch^~^[p^'yxY^'^~^{PylyY^'^l"'^ where ifji, (1 — il)^) are supported in {|x — y| > 7} and 
— y| < 27} and 7 > hp~^ . 
Furthermore, since 

(2.62) \Q,EQy\ < Ch'-'p'J^-'r.'p'J^ 

due to the standard arguments, contribution of the pair (/ — il)x)Qx and Qy to the error 
does not exceed Ch^~^^ p^~^ Pyj~^'yyj'^~'^ Plugging 7 = hp~^ we estimate the contribution 
of the pair Q^, Qy by 

(2.63) Ch'-'-{p,l.Y^'~\PylyY^'p: + Ch'~'-p~'^''i;'p'yYy 

which is larger than (2.61). 
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In these estimates we do not need non- degeneracy condition and therefore as (y, rj) and 
(x, ^) are given we can take 



(2.64) p^ = py = |x-y|'^ + |^-ry|,7^ = 7y|x-y| + |^-7/ 



l/a 



where p = 7*^ on the corresponding step of our analysis. Then as (z, Q) are fixed contribution 
of {|x — z| < 7, |y — z| < 7, 1^ — CI < P< I'? ~ CI ^ P< k ~ y| + I'C ~ ^ ^7} to the error 
does not exceed this expression 

(2.65) Ch^~'^-''{p-iY-^p^ + Ch^-''-^[p-iY-^p'' 

where the second term is less than the first one. 

Then the total contribution of the zone in question to the error does not exceed 



(2.66) C/?^-^ / / / 7'^«-'^-i c/yc/r/ 7-M7 




where equation is taken over {7 > 7x} and the integral in question is equivalent to Mh^ 
where M = 1 as a{n — 1) > 1, 



(2.67) jjj 

as a{K — 1) = 1 and due to (2.47)n M x 1 as well, 

(2.68) jjj ^^"'^""^ ^^^^ 

as (t(k — 1) < 1, and on each step of the induction we already proved that M x 1. □ 
2.6 General Microhyperbolic Case. II 

Let us consider matrix operator. Let Aj(x, ^) be eigenvalues of its principal part. Then 
iVx.^Ajl < c and microhyperbolicity with respect to I means that 

(2.69) |A,(x,OI<eo =^ > eo Wj. 
Let us consider zone 

(2.70) ^■ = {(x,0:|Ay|<min|A,|} 
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and let us define here 

(2.71) 7 =^ min \\k\ + 
and p = 'J- Consider zone 

(2.72) {7> |x-y| + |e-r/|+7} 

and let us rescale x x/'j, ,^/7, Aa- ^ ^k/l, h ^ /'/7^ preserving microhyperbolicity 
condition (1.15) and simultaneously making operator with > 1 for /c 7^ _/ and therefore 
analysis of this operator is not different from the scalar one. Unfortunately we cannot 
use non-degeneracy conditions of subsections 2.4-2.5 which would not survive this, but 
microhyperbolicity condition survives and we assume that {1.15) is fulfilled. 

Then as the main part of the asymptotics is given by the standard Weyl expression 
(2.27) — (2.29), the contribution of zone (2.72) (intersected with {7 > C07}) to the remainder 
does not exceed 



(2.73) Rj = [ c{hr'y~'~^r-^' dpj X ch'-'~- [ r'-"^ dpj 

.Jzjn{'r>C^} -'5:,n{7>7} 

with Tj = {(x, ^) : \j = 0} and dp; = dxd^ : dXj density on it. 

Let us fix 7 = Ch^^^. Then in the complementary zone U/(^j{|Aj| + |A/<| < C7} one needs 
just to make a rescaling xi-^x/7, ,^1-^,^/7 which sends /? to 1 and no microhyperbolicity 
condition would be needed and the contribution of this zone would not exceed 

(2.74) R'j;^ = C/)-^-"/2 mes{|Aj| + |A/,| < Ch^^^}. 

So, the total contribution of zone UjUj to the remainder is given by J^j + J2j k j^k ^jk- 
Assuming that 

(2.75) pj{Tj : |A/,| < t) + mes{\Xj\ + |A/,| < t} = 0{f). 

we get under additional assumption r + k, > 2 (which is always fulfilled as r > 2) that 
Rj = 0{h^-''~^) while RL = 0(/7^) with 



(2.76) 




which is 0{h^ as well. 
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On the other hand, as r + k < 2 we get that Rj = (/?''), Rji^ = 0{h^) with q given by 
(2.77). 

Finally, as r + k = 2 we get Rj = 0{h^-''~^\ log h\) and Rjk = 0{h^-''~''). 

Assume temporarily that no more than two eigenvalues can be close to simultaneously. 
Then we are already done since in the zone complimentary to (2.72) we redefine 7 = e(|x — 
y| + — r/l) and apply the same rescaling as before and one does not need microhyperbolicity 
condition. 

Let us apply induction by m assuming that no more than m eigenvalues can be close to 
simultaneously. Then we can define on each step 

(2.77) 7(x,0 = e max min |A/,(x, ^1 + 7 

J:#J=m k<^J 

and repeat all above arguments. We will arrive to 

Theorem 2.20. Let conditions [0.9), {0.7), {1.15) and {2.75) he fulfilled. Then the standard 
Weyl asymptotics {2.27) — {2.29) holds with the remainder estimate 

(i) which IS 0(/)i-^-('"-i)'^) as r + K>2; 

(ii) which is 0{h'^) with q defined by {2.76) as r + qn < 2 and 0{h^^'^~'^\ \ogh\) as r + n = 2. 

Remark 2.21. Condition (2.75) is fulfilled provided Aji^ = {Xj = = 0} are smooth 
manifolds of codimension r and \Xj\ x \Xk\ x dist((x, C,), ^jk) in its vicinity; this assumption 
should be fulfilled for all j ^ k. 
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